Spin density wave in oxypnictide superconductors in a three-band model 
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The spin density wave and its temperature dependence in oxypnictide are studied in a three- 
band model. The spin susceptibilities with various interactions are calculated in the random phase 
approximation(PPA). It is found that the spin susceptibility peaks around the M point show a spin 
density wave(SDW) with momentum (0, 7r) and a clear stripe-like spin configuration. The intra- 
band Coulomb repulsion enhances remarkably the SDW but the Hund's coupling weakens it. It is 
shown that a new resonance appears at higher temperatures at the P point indicating the formation 
of a paramagnetic phase. There is a clear transition from the SDW phase to the paramagnetic 
phase. 

PACS numbers: 74.20.-z, 74.25. Jb, 74.72.-h 



I. INTRODUCTION 

The high temperature superconductivity in the newly 
discovered oxypnictides, LnFeAsO (Ln = La,Pr,Ce,Sm 
), has attracted great attention aiming to identify the 
mechanism of superconductivity in these materials [l|. 
Recently the transition temperature T c is dramatically 
raised from 26A' to AZK[2). In addition to the high T c , 
these materials display many other interesting proper- 
ties. The most interesting phenomena is the presence 
of competition between the magnetically ordered ground 
statesQ of spin density wave(SDW) and superconduc- 
tivity, but a controversy model explains the supercon- 
ductivity by means of antiferromagnetic spin fluctuation 
in LaFeAsO[4]. Pure oxypnictide is not superconducting 
but shows an anomaly at about 150K in both resistivity 
and dc magnetic susceptibility P, Q. Both experimen- 
tal and theoretical evidences show that the anomaly is 
caused by the SDW inst ability d H 0|. 

It is shown by the first principle calculations that the 
band structure of LnFeAsO near the Fermi surface (FS) 
is formed by a hole-like pocket centered around the T 
point and an electron-like pocket around the M point in 
the extended Brillouin zone(BZ)(one Fe atom per unit 
cell) [||-[l(| . A strong FS nesting effect exists between 
the hole and electron pockets with commensurate wave 
vectors, (it, 0) and its symmetric ones. This leads to a 
strong SDW instability, and is believed to cause the the 
' anomaly at 150K. 

Raghu et al calculated the SDW in the iron oxypnic- 
tides within a minimal two-band model in the phase ran- 
dom approximation(RPA). They found that the SDW 
is enhanced significantly by the intra-band Coulomb re- 
pulsion. The influence of the inter-band Coulomb inter- 
action and the Hund's coupling, however, have not yet 
been fully studied in literatures. Different researchers se- 
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FIG. 1: (a) A schematic diagram showing the hopping param- 
eters between d xz (d yz ) and d xy . (b) The hopping parameters 
between the nearest and next nearest d xy orbitals. 

lcct different groups of interaction parameters, but the 
relation between them has not been revealed. In addi- 
tion, the temperature dependence of the SDW has hardly 
been theoretically considered. As pointed out by Lee and 
WenfTTj. a three-band model reproduces more accurately 
the band structure near the Fermi surface of oxypnic- 
tides. In this work we study the SDW in a three-band 
model in the RPA and the temperature dependence of 
the spin susceptibility. The magnetic instability is stud- 
ied for a wide range of interaction parameters. We found 
a new resonance of the spin susceptibility around the T 
point at higher temperatures, indicating the formation of 
a paramagnetic phase. There appears a transition from 
the SDW phase to the paramagnetic phase when temper- 
ature increases. 



II. MODEL HAMILTONIAN 

The FeAs layer in LaFeAsO forms a square lattice, 
where Fe ions locate on the lattice sites and an As ion 
sits at the center of each square. Various band structure 
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calculations showed that the main contribution to the 
density of states near the FS comes from d xz , d yz and 
d xy of the five 3c? orbitals of Fe atoms [ll|. The left two 
orbitals are far apart from the FS. 

In the three band model, the hopping terms between 
and d xy as illustrated in Fi g.1 are included in 
the Hamiltonian. They are written as [ill ll2T ] 

ff0=5>L^k*fca, (1) 

where the three-component field is defined as ^>ka = 

(d xzcr (k),d yzcT ('k),d xy<T ('k)) T and the Matrix M k is given 
by 

(£n(fe) £12 (fe) £13 (fe) \ 

£2i(fe) e 22 (fe) e 23 (fc) , (2) 

£31 (fe) £32(fe) £33(k) J 

with elements 

£11 = — 2tj cos k x — 2t 2 cos k y — 4i 3 cos k x cos k y 
£22 — ~2t 2 cos k x — 2ti cos k y — At 3 cos k x cos k y 
£33 = — 2tj cos k x — 2tj cos fey — 4ig cos k x cos fe y 
ei2 = £21 = — 4^4 sin fca; sin fc y 
£13 = -£31 = -2^5 sin fez 
£23 = -£32 = -2it 6 shxk y . 

The hopping parameters are set to t\ — — 1.0 (~ 
0AeV),t 2 = 0.7, t 3 = -0.80, i 4 = 0.6, i 5 = = 
-0.35, i 7 = -0.3, i 8 = 0.2, fx = 1.15, in units of [V2}. 

The energy dispersion is plotted in Fig. 2(a) in the 
extended BZ and 2(b) along the T X M ->■ T in 
the folded BZ(two Fe atoms per unit cell). Fig. 2(d) plots 
the density of states of the band structure. It is seen that 
there are two dominating Van Hove singularities near the 
Fermi level, close to the hole and electron pockets. The 
density of states is significantly different from the two- 
band model, of which the Van Hove singularities are more 
distant from the Fermi level. A stronger nesting effect is 
expected in this three-band model. The third band is 
mixed with the two conventional bands and thus should 
have great contribution to the spin susceptibility. 

III. THE SPIN SUSCEPTIBILITY 

The static spin susceptibility in the non-interacting 
case is given by x (0) (q) = E«' X?z'(q) with 

xk(q) - m E fc /( trj ( -^ q ' ,} I < * + g, > l 2 (3) 

where I, I' = 1, 2, 3 are band indexes, /(e) = \/{e^ e -^ + 
1) is the Fermi distribution function, /3 = 1/feT, and Em 
and |k, I > are the eigenvalues and eigenvectors respec- 
tively of the Hamiltonian matrix @. 

We fix kT = 0.02(- 93K) for a finite lattice with 64 x 
64 k meshes in the extended BZ to calculate the static 
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FIG. 2: (a) Energy dispersion in the unfolded BZ of the three- 
orbital model with U = -1.0, t 2 = 0.7, t 3 = -0.80, U = 
0.6, ig = U = -0.35, tr = -0.3, t& = 0.2,^ = 1.15 along the 
path (0,0) — s> (7r,0) ->■ (71-, 7r) — > (0,0). (b) Energy dispersion 
in the folded BZ along r = q = (0, 0), X = q = (§, §),M = 
q= (71", 0). (c) Energy dispersion on the (k x ,k y ) plane, (d) 
density of states of the three-band model, where the dashed 
line represents the Fermi level. 



spin susceptibility, which is shown in Fig.3(upper). The 
static spin susceptibility shows great peaks at (0, ±7r) 
and (±7r, 0). This indicates that a collinearly-striped an- 
tiferromagnetic (AFM) order phase, a spin density wave 
(SDW), exhibits in oxypnictide. This feature is in agree- 
ment with the neutron scattering measurements [M 0TLl3 ■ 
The SDW comes from the strong nesting between the 
hole pocket at the T point and the electron pocket at the 
M point in the extended BZ. It is observed experimen- 
tally that this SDW peak is significantly suppressed by 
F doping Q. This is reasonable because a down-shift of 
the Fermi level tends to reduce the size of electron pocket 
and enlarge the hole pocket thus to suppress the nesting 
between them. 

Now we consider the interactions: the intra-band 
Coulomb repulsion U, the inter-band coulomb interaction 
U', the Hund's coupling J. The interaction Hamiltonian 
is written as 

Hmt = U ^^nufnui + U' nungi 

il i,l^V 

+ JJ2 S U- S '"' ( 4 ) 
i,ljtl' 

In the RPA[l4l|. the spin susceptibility, x s (q)i an d the 
charge-orbital susceptibility, x c (q), with interactions are 
given by 

X s (q) =[/-f/V(q)]-V(q), (5) 
X c (q) =[i + [/ c X°(q)rV(q), (6) 



FIG. 3: The static spin susceptibility x (o) versus q(upper). 
The static spin susceptibility x(q) i n the RPA(lower) with 
U — 3.0, J = 0., kT — 0.02. 



FIG. 4: The determinants of magnetic instabilities for J vs. U 
and U' vs. 17. The contour line —0 — 0— gives the parameter 
groups of instability. 



where X°( ( l) is a 3 x 3 matrix with elements denned in 
Q and II s ( c ) are the interaction matrices 



II s = 



U -J -J 
-J U -J 
-J -J u 



,u c = 



U 2U' 2U' 
2U' U 2U' 
2U' 2U' U 



It is seen from ([5J and (J6j) there will appears a magnetic 
instability when the following relations are satisfied: 



£fei[JTE/ s(c) X°(q)] =0 



(7) 



To show the magnetic instability at the M point the de- 
terminants (|7|) at the SDW momentum, (0, 7r), are calcu- 
lated for different interaction parameters and are plotted 
in Fig. 4. The parameters on the contour lines labeled 
by "— 0— " in these two diagrams give zero determinants 
thus lead to magnetic instabilities. For a fixed U the spin 
susceptibility decreases with increasing Hund's coupling 
J thus reduces the SDW peak at the M point. This re- 
sult is quite different from that made by Raghu et q?(l5| , 
who found stronger spin fluctuations with increasing J. 
The lower diagram in Fig. 3 shows the spin susceptibility 
close to the magnetic instability with U = 3.0, J = 0.0, 
where the SDW peak is much enhanced due to the in- 
traband Coulomb repulsion relative to that of the bare 
spin susceptibility. When J increases, however, the SDW 
peak drops significantly until the SDW phase fully dis- 
appears. It is interesting to notice that, when U < 5, 
a ferromagnetic coupling (J < 0) is beneficial to the 
formation of the SDW on the M point. On the other 
hand, the charge-orbital instability appears in the region 
—4.2 < U' < 0. That is, when a Coulomb attraction 



exists between different bands a charge-orbital instabil- 
ity may occur. Apart from this region the charge-orbital 
susceptibility also drops rapidly. It is worthy to notice 
that the charge-orbital susceptibility depends weakly on 
the value of U, the intra-band repulsion, except for the 
case of very strong inter-band repulsion. 

The SDW phase has striking temperature dependence. 
It is found that at higher temperatures a new resonance 
appears at the T point, which starts to increase rapidly 
at a temperature kT ~ 0.08, as shown in Fig. 5. This 
resonance corresponds to the formation of a paramag- 
netic phase in the material. The SDW peak at the M 
point drops to a nearly stable value at this tempera- 
ture. This result is qualitatively in agreement with the 
experimental observation for the stripe-like AFM phase 
in LaOFeAs, which forms under a temperature 134K[l6j. 
At kT c = 0.12 the intensity of paramagnetic phase starts 
to surpass that of the striped AFM phase. It is found 
that temperature changes hardly the bare spin suscepti- 
bility but smooths its distribution in the BZ. Reducing 
the intra-band coupling parameter U, the relative inten- 
sity between the paramagnetic phase and the SDW phase 
reduces significantly. Hence it is confirmed that this new 
resonance comes from the intra-band coupling. There- 
fore, this temperature dependence of the spin suscepti- 
bility provides detailed information for the band struc- 
ture of the oxypnictide material. More experimental ev- 
idences for this temperature dependence are expected. 
At very low temperatures the spin susceptibility showed 
more abundant structures but requires further stringent 
calculations. 

We tried to solve the superconducting gap equation 
by including a pairing coupling term into the interac- 
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FIG. 5: Spin susceptibility with temperature T = 0.2(upper) 
and temperature dependence of the intensities of spin sus- 
ceptibility at the F point (circles) and M point(stars) in the 
BZ(lower). Both results take U = 3.0, J = 0, V = 1.0. Lines 
are guides for the eye. 

tion Hamiltonian but failed to find a stable gap func- 
tion with some symmetries, such as s, s~, d x 2_ y 2, d xy , 



etc. This failure also exists in some literatures, e.g., Yao 
et al found only eigenvalues 0.1-0.4 with the d xy — wave 
symmetry [13] ■ But Yanagi et al claimed that some 
groups of interaction parameters give eigenvalue unity 
for this gap symmetry in the same material (l8|. There- 
fore, the gap symmetry remains controversy and requires 
further studies. 

In conclusion, a three-band model is set up to repro- 
duce the band structure near the fermi surfaces of oxyp- 
nictidc. It shows a hole pocket around the T point and 
a electron pocket around M point in the extended BZ. 
The spin susceptibility with various interactions are cal- 
culated in the random phase approximation. It is found 
that the spin susceptibility peaks around the M point, 
showing a SDW with momentum (0, tt) and a clear stripe- 
like spin arrangement. The intra-band coupling enhances 
remarkably the SDW but the Hund's coupling weakens 
it. The interaction parameters are determined for the 
magnetic instability in oxypnictide. Finally the temper- 
ature dependence of the spin susceptibility is studied. It 
is found that a new resonance appears at higher tem- 
peratures at the r point indicating the formation of a 
paramagnetic phase. There is a transition between the 
SDW phase and the paramagnetic phase. 
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